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Abstract 

We study Stokes phenomena of the k x k isomonodromy systems with an arbi- 
trary Poincare index r, especially which correspond to the fractional-superstring (or 
parafermionic-string) multi-critical points (p, q) = (1, r — 1) in the fc-cut two-matrix 
models. Investigation of this system is important for the purpose of figuring out 
the non-critical version of M theory which was proposed to be the strong-coupling 
dual of fractional superstring theory as a two-matrix model with an infinite number 
of cuts. Surprisingly the multi-cut boundary-condition recursion equations have a 
universal form among the various multi-cut critical points, and this enables us to 
show explicit solutions of Stokes multipliers in quite wide classes of {k, r). Although 
these critical points almost break the intrinsic symmetry of the multi-cut two- 
matrix models, this feature makes manifest a connection between the multi-cut 
boundary-condition recursion equations and the structures of quantum integrable 
systems. In particular, it is uncovered that the Stokes multipliers satisfy multiple 
Hirota equations (i.e. multiple T-systems). Therefore our result provides a large 
extension of the ODE/IM correspondence to the general isomonodromy ODE sys- 
tems endowed with the multi-cut boundary conditions. We also comment about 
a possibility that N = 2 QFT of Cecotti-Vafa would be "topological series" in 
non-critical M theory equipped with a single quantum integrability. 



*ctchan@thu.edu.tw 

tirie@phys.cts.nthu.edu.tw 

*d95222008@ntu.edu.tw 



Contents 



1 Introduction [ll 

2 Summary of fractional-superstring critical points 5 

2.1 Definition of tlie multi-cut matrix models 5 

2.2 The ODE system and the weak-coupling results § 

2.2.1 Two kinds of critical points and their spectral curves 6 

2.2.2 Superposition of perturbative string theories 

2.2.3 The multi-cut geometry and Stokes phenomena 

2.3 Stokes phenomena in the multi-cut critical points 

2.3.1 Asymptotic expansions 

2.3.2 Stokes sectors 

2.3.3 Profiles and Stokes multipliers 

2.3.4 Formula of profile components 

2.3.5 Three basic constraints 

3 The multi-cut boundary conditions (BC) 

3.1 The multi-cut BC recursion equations 

3.1.1 The boundary conditions and the shift rule 

3.1.2 The theta parametrization 

3.1.3 The multi-cut BC recursion equations 

3.2 Relation to discrete Hirota dynamics 

3.2.1 Hirota dynamics for non-critical string/M theory 

3.3 Explicit solutions of Stokes multipliers 

3.3.1 The case of k = rm + 1 

3.3.2 The case of A; = rm + I, {1 < I < r) 

3.3.3 The case of /c = r(m -|- 1) 

4 Conclusion and discussions 
A Sketch of a proof for Theorem [1] 



1 Introduction 

Non-perturbative aspects beyond perturbative string theory still remain some of the 
missing pieces in our current understanding of string theory. In early study in the '90s, 
the following question about non-perturbative string theory was investigated: what is 
the form of non-perturbative corrections to the perturbative expansions of string theory? 
Among various investigations, non-critical string theory [Ij played an important role since 
it is non-perturbatively formulated by solvable matrix models [2l-H^ and this fact makes 
it possible to see higher-order behavior of the string-theory free energy in strong string 
coupling constant gstr oo [T5l[l6l|211|22l|2ll|25]. In particular, it was found that non- 
perturbative corrections to the perturbative amplitudes are given by order of (9(e~*/^""'), 
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i.e. open-string degree of freedom 



oo 
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perturb, part non-perturb. corrections 

and this fact lead to a key idea of D-branes [15]. The identification of D-branes was then 
confirmed in comparison between the Liouville-theory calculation [T]H6] - [55] . i.e. ZZ-branes 
|51j . and succeeding calculations of matrix models [3TI435ll37ll39] . On the other hand, since 
they describe non-perturbative corrections in the perturbative asymptotic expansion, the 
relative weight factors {^/j^, called D-instanton chemical potentials (fugacities), remain 
to be ambiguous parameters in perturbative string theory |21j . 

A basic origin of the ambiguity is the fact that these parameters are subjected to 
discrete corrections in various analytic continuations of string-theory backgrounds and 
that the formal expansions of analytic equations can be performed without specifying 
the direction from which we expand the series (See Figured]). From the viewpoint of 
physics, this kind of phenomena make it possible to describe various different perturba- 
tive backgrounds in a single analytic theory. Mathematically these behaviors are known 
as Stokes phenomena, a basic analytic property in the asymptotic analysis. Determining 
the Stokes phenomena for the asymptotic expansions in the string coupling constant, 
Eq. fll.ll) . (i.e. how the fugacities {Oi}j jump in analytic continuations) can be per- 
formed by resurgent analysis of the perturbative asymptotic expansion Eq. (II. ip . which 
was analyzed recently |56l|57]. On the other hand, the non-perturbative ambiguities are 
information completely out of the asymptotic expansions, and therefore we have to come 
back to the original definition of matrix models and need to see additional principle in 
order to determine the D-instanton fugacities in string theory. This is the question we 
investigate in this paper: what is the principle to choose physical D-instanton fugacities? 
A reason why we consider these non-perturbative problems so significant is because we 
believe that it is this principle which binds various string-theory vacua together to for- 
mulate non-perturbatively completed string theory and consequently which guarantees 
the string dualities in a consistent way. Therefore, this would be the direction in which 
we can quantitatively tackle the string-theory landscape from the first principle. 

In study of this issue about the non-perturbative ambiguity, non-critical string theory 
is again expected to play an important rule. In non-critical string theory, the non- 
perturbative ambiguities are identified as integration constants of string equations |21j . 
At the first sight, therefore it may seem plausible to say that the D-instanton fugacities 
are arbitrary free parameters in the non-perturbative theory. However, it is also true that 
there are several systems which prefer "physical values" rather than arbitrary values [58] , 
and these physical values are implied in the matrix models. This kind of study has been 
investigated from various viewpoints [58]ff2] . In this paper, succeeding the previous study 
of the authors [72] , we investigate this issue in the framework of the multi-cut two- matrix 
models (SUES] [131I111E2]. 

The multi-cut two-matrix models are the two-matrix version of the multi-cut one- 
matrix models [73] and realize the multi-cut critical points which naturally extend the 
two-cut critical points of the matrix models [751480] . In particular, these critical points 
are controlled by mult i- component KP hierarchy [52] and this fact makes possible var- 
ious quantitative analysis of the matrix models [15 | 1 ^ 172]. What is more, the two-cut 
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Figure 1: The asymptotic expansions can be different up to the D-instanton fugacities {^/} ^ depending 
on from which direction of the string coupHng gstr (— ^ X e*'') we expand tlic free energy (or partition 
function). That is, in this figure, there are eight different sets of the D-instanton fugacities {d'f ''} j (a — 
1, 2, • • • ,8) depending on the angular domains of the string couphng gstr- This mathematical property 
is called Stokes phenomenon. Therefore, determining Stokes phenomena means determining how the 
fugacities jump in analytic continuations crossing the domain wall (e.g. O^f'' — 0\^^ in crossing the 
positive real axes). Depending on the relative magnitude of perturbative free energy T^^'^igstv) and the 
D-instanton fugacities, various perturbative string theories emerge as a dominant saddle point. This 
kind of phenomena can occur in any analytic continuations of the string-theory moduli spaces. 



critical points were found to describe minimal type superstring theory [8T1 - I83] . and, as 
an extension of this, a special kind of the multi-cut critical points (which is also shown 
in Section I2.2.ip were proposed to describe minimal fractional superstring theory (or 
minimal parafermionic string theory) [73j. An interesting new feature found by quanti- 
tative analyses of these critical points is the fact that these multi-cut systems generally 
include a number of perturbative vacua in its string-theory landscape |13] and also in- 
clude various perturbative string-theory sectors in its wavefunctions [H]. In particular, 
the fractional-superstring critical points are found to describe a superposition of various 
minimal fractional superstring theories, and the chain of string theories is then inter- 
preted as the extra dimension of the non-critical M theory [H]. This means that, in the 
case of an infinite number of cuts, there naturally appears a three-dimensional universal 
strong-coupling dual theory which realizes the philosophy proposed by Hofava-Keeler for 
the non-critical M theory In this way, various new non-perturbative phenomena of 
string theory were revealed in the multi-cut critical points, and importantly these phe- 
nomena are fully described only after we determine the non-perturbative ambiguities, 
i.e. they are governed by the information out of the perturbative description of string 
theory. Therefore, it is clear that the multi-cut matrix models provide a fruitful field 
in studying this issue, and the investigation of this system is expected to answer to the 
following question: what is the non-perturbative completion of string theory? 

As a natural framework to study the ambiguity, we consider the isomonodromy de- 
scription [SSIIHS] [HUSH] [SZ] for the multi-cut critical points. Mathematically, the system 
is a. kxk isomonodromy system with an irregular singularity of Poincare index r (= q + 1) 
which corresponds to the fc-cut multi-critical points {p,q) = (1, q) in the multi-cut ma- 
trix models. In the isomonodromy study of matrix models, the D-instanton chemical 
potentials are mapped to Stokes multipliers of the linear systems in the spectral parame- 
ter ( (i.e. the Baker-Akhiezer systems Eq. fl2.14p ) through the Riemann-Hilbert procedure 
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(See ^87] ). Roughly speaking, the Riemann-Hilbert procedure provides an integral ex- 
pression of the string free energy J^((7str) which is given by the Stokes multipliers of the 
ODE systems. For instance in the case of the 2-cut {p,q) = (1,2) critical point, the RH 
integral is give by the Stokes multiphers {sn,i,j}n,i,j of the corresponding ODE system as 



where t is a parameter of the string-theory moduli space (i.e. the worldsheet cosmological 
constant) 111 In this way, the information of the non-perturbative ambiguities is translated 
into the Stokes data of ordinary differential equation (ODE) systems as the information 
which does not depend on the direction of asymptotic expansion (shown in Figure 
Although the Stokes phenomena in these ODE systems are relatively simpler than those 
of non-linear differential equation systems like Painleve equations, less is known in the 
case of large values of k and r compared with the isomonodromy systems of the second 
order Painleve series |87j . 

In the previous study of the authors [72j , the physical constraints (required from the 
matrix models) on the Stokes multipliers are formulated in this isomonodromy framework. 
As is reviewed in Section I2.2.3[ the requirement is given as the boundary conditions for 
physical cuts in the spectral curves. Interestingly these physical constraints enable us to 
show explicit solutions of Stokes multipliers which satisfy high-degree algebraic equations 
called monodromy free (or cyclic) conditions (Eq. f l2.65p ). which are generally quite hard 
to solve without any help of the physical boundary conditions. 

In this paper, on the other hand, we extend this boundary condition to wider classes 
of the multi-cut critical points (general number of cuts. A;, and general Poincare index, 
r), especially the critical points which correspond to minimal fractional superstring (or 
parafermionic string) theory |73j. A major difference from the critical points of the pre- 
vious study is lack of the Z^-symmetry which is intrinsic in the multi-cut two-matrix 
models. Therefore, the number of independent Stokes multipliers are k times more than 
that of Zfc symmetric critical points, and this becomes a huge predicament to direct anal- 
ysis of the Stokes phenomena. However, we will show that by this feature it becomes 
manifest that the multi-cut boundary conditions are equivalent to discrete Hirota equa- 
tions (i.e. T-systems) of quantum integrable systems. That is, the physical section of 
D-instanton chemical potentials are expressed by quantum integrability which would be 
an integrable structure of strong-coupling dual description, i.e. non-critical M theory. 

Organization of this paper is as follows: In Section [2|, basic facts about the fractional- 
superstring critical points are summarized. In Section [31 the multi-cut boundary condi- 
tions are studied. In particular, the T-systems of the Stokes multipliers are pointed out 
in Section [3l2l Some explicit solutions are also shown in Section [3731 Section [His devoted 
to conclusion and discussion. In Appendix [Al a brief proof of Theorem [H is presented. 



^One can find an explanation of this integral in [72] and |87j . 

^For example in |87| . by evaluating this integral, they show explicit expressions of how the figacities 
{Oi} J (i.e. for each angular domain (a = 1,2, • • •) in Figure [1]) can be expressed by the Stokes 

data {sn.i,j}n.i,j of the ODE system of the Painleve equations. 
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2 Summary of fractional-superstring critical points 



In this section, we summarize basic facts about the fractional-superstring critical points 
[75], especially focusing on the results from the matrix- model analysis I131I111E2] , includ- 
ing the Stokes phenomena of the multi-cut matrix models which are generally developed 
in 



2.1 Definition of the multi-cut matrix models 

The multi-cut two-matrix models are given by two-matrix models, 

-N tr [Vi iX)+V2 iY)-XY] 



Z= / dXdYe-''''r'^''>+''^'>-^'i, (2.1) 

iMjv(C('=))xMjv(C(*)) 

with integration over the following special N x N normal matrices, M^^C^^'^): 

M,v(C('')) = {f/diagti(A,)?7T|?7 G U{N), \i e C^''^ C C}, (2.2) 

where the contour of the eigenvalues, C^''\ are defined as 

fe-i 

C^'^) = IJ u;"M, a; = e'"*i, (2.3) 

n=0 

in the case of the fc-cut two-matrix models. 

An important ingredient of this system is the bi-orthonormal polynomial systems [55] 
defined by the following inner product, 

rfxf/i/e-^[^^(^)+^^(^)-^^] an{x) f3m{y) = 5„,„., (2.4) 

between the polynomials, 

a„ (x) = -1= (x" + ■ ■ ■ ) , p^(y) = -}-(yn + ...y (2.5) 

Here {hn}nez+ are normalization coefficients. An interesting fact about the ortho normal- 
polynomial system is the following expression by the matrix integral [llj: 

^n{x) = ^(det(x-X)>_, /3^{x) = ^(det(y-r)>^^^, (2.6) 

where the expectation value ■ • is given by n x n truncated matrices Mn{C^^^): 

{0{X,Y)),,, = [ f/Xf/ye-^*4^iW+^2(^)-^^] 0(X,F). (2.7) 

Another important ingredient of this system is the resolvent operator R{x) which 
contains the spectral information of the matrix integral: 

i?(a;) = l/tr^-\ ^ / dX^^ (N-^oo). (2.8) 
^ ^ N\ x-X / Jew x-\ ^ ^ ^ ^ 
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In particular, the support of the eigenvalue density function p(x) corresponds to discon- 
tinuity of the resolvent function: 

7iip{\) = R{\ + ie) - R{\-ie), A G supp (p) C C^^^ (2.9) 

Since the orthonormal polynomials are also written with the resolvent operator: 



anix) 



N 



det(x -X)) 



: exp 



N / dx'R{x) (n~ X), (2.10) 



the above discontinuity should be observed also in the orthonormal polynomial of the 
large N limit with n/N — )■ 1. 



2.2 The ODE system and the weak-couphng results 



We move on to the continuum {p, q) minimal fractional superstring theory [73] by taking 
the scaling limit of the bi-orthonormal polynomials. The scaling limit is given by a lattice 
spacing a as follows: 



M(t;C)\ 



^orth(t;C) = 



A{t-X) = aiu;^/2[(-l)"a,„+,_i(x)e-^^^(^) 



where the scaling valuables are introduced as 



(2.11) 



x = a^uj^'^(, N=gstra^ 0. 

The detail of these scaling limits is discussed in 
a parameter e (= 0, 1), which are related to the two choices of fractional-superstring {p, q) 
critical points. Consequently, the recursion equation of the orthonormal polynomial: 



kn , p+g-i \ 

— =exp(-ta 2 j ^ 1, 

a->0). (2.12) 
As will be mentioned later, there is 



n ar, 



(X 



d 

N-^—an{x) 
ox 



(x) 



(2.13) 



becomes the following k x k differential equation system [39] : 

Cvl>(t; C) = P(t; d) vl/(t; C), ^?str^^(t; C) = Q{t; d) vl/(t; C). 



(2.14) 



which is governed by /c-component KP hierarchy |89]. Here we note that the scaling or- 
thonormal polynomial \&orth(^i C) of Eq. (12. lip is a special solution to the Baker- Akhiezer 
function system Eqs. f l2.14p . 



2.2.1 Two kinds of critical points and their spectral curves 

The critical points of the multi-cut two-matrix models are classified by these lax operators 
(P, Q). Fractional-superstring (p, q) critical points are originally discussed in [73] as the 
critical points which realize the operator contents of {p, q) minimal fractional superstring 
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theory. From the quantitative analysis of these critical points [43] (which is the multi-cut 
extension of |23|), it was found that there are two kinds of fractional-superstring critical 
points which are called co^^'^ -rotated models and real- potential models. The parameter e 
in Eq. ( ]2.12p then is related with these two choices of the fractional-superstring critical 
points: 



: cj^/^-rotated models 

1 : real-potential models 



(2.15) 



These two kinds of solutions are perturbatively distinct when k G 4Z; non-perturbatively 
always distinct. Here we summarize the results of critical-potential analysis |43j and the 
weak-coupling spectral curves 



1. The first kind is derived from u^^'^ -rotated-potential models [43j, and given as 

p q 

p{t- d)=Td^ + Y^ i/f {t) d^-"", Q{t- 9) = r 9^ + ^ iyf {t) 9^-' 



n=l 



n=l 



(2.16) 



with the shift matrix F, 



/O 1 




1 



(2.17) 



0/ 



and the real functions H^^^\t) and H^'^^\t). The weak coupling analysis gstv 
results in the following simultaneous eigenvalues of the Lax operators P{t; d) ^ 
die.g';^,{p(^\t;z)),Qit;d)^dmg';^,{Q(^\t;z)) [44]: 



Cosh solution: 

pO) = APcosh(j3r + 27Ti^-^), 

k 

Sinh solution: 

P(^) = APsinh(pT + 2m^^-^), 

K 



qO) = A«cosh(gr + 27ri^-— ); (2.18) 

k 



Q(^) = A^sinh(gr + 2 



.J 



m- 



-), (2.19) 



with 

\ = tT^\ 2 = A-(P+^)^,t,t9t = cosh(r). 

Their algebraic equations are expressed a^ 

Cosh solution: F(P, Q) = Tp{Q/X'^) - T^^PjX^) = 0; 

Sinh solution: F(P, Q) = Tp{-iQ/\^) - (-l)^Tg(-iP/A^) = 0. 



(2.20) 

(2.21) 
(2.22) 



^rp(cosh(T)) = cosh(pT) stands for the Chebyshev polynomial of the first kind 
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Here the indices {p, q) are identified with the labehng of minimal fractional super- 
conformal models [73] and defined by 

{p,q) = {kp,kq), (2.23) 
with introducing the following integ ers dq—p, k and Tj^ 

dg-p = g.c.d.{q-p,k), k = kxdg-p, q - p = i] x dg-p. (2.24) 

2. The other kind is derived from real-potential models [13] and given by replacing the 
shift matrix F by the twisted shift matrix r(i.eai) as follows: 



P{t; d) = r(real) d^ + Yl H^r^'-'^^t) Qit'^ 9) = r(,eal) 0^ + Yl '"^H^) 9' 

n=l n=l 

(2.25) 

with the twisted shift matrix r(i.eai), 

\ 



,q-n 



/O 1 

1 



(real) 



-1 



1 



(2.26) 



and the real functions H^''^\t) and H^^^\t). The weak couphng analysis gstr 
results in the following eigenvalues of the Lax operators P(t; d) ~ diag^^^ (P*^-'^(t; z)^ , 
Q(t;9)^diag^ti(Q(^)(t;z)) [MJ: 



Cosh solution: 

p(j) = A^cosh(pr + 2 
Sinh solution: 

P(^) = A^sinh(pr + 2 



m- 



TXl- 



.2 j - 1 

2k 

..2 j-l 
2k 



), g(i) = A«"cosh(gr + 27ri^:^); (2.27) 



) , gO) = sinh (qr + 2m-^—^) , (2.28) 
' 2k 



with the parameters in Eqs. f l2.20p . Their algebraic equations are expressed as 

Cosh solution: F(P, Q) = Tp{Q/X^) - (-1)^ Tg{P/\^) = 0; (2.29) 
Sinh solution: F{P, Q) = Tp{-iQ/\^) - (-l)^^^+^T,(-zP/A^) = 0. (2.30) 

Here rj is defined in Eq. f l2.24p . 
2.2.2 Superposition of perturbative string theories 

A nontrivial feature of these explicit expressions [S] is that these algebraic equations 
F(P, g) = are generally factorized into several irreducible pieces of the algebraic curves 
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FW(P,Q) = 0, (j = l,2,---,A;)g 



LlJ+i 

F{P,Q)= n F^'\P,Q) = 0, (2.31) 

and these irreducible pieces are equivalent to the algebraic curves of the fractional- 
superstring critical points (of the matrix models) with the smaller number of cuts. For 
instance, some are equivalent to the curves of minimal bosonic string theory and some 
are those of minimal type superstring theory. 

The authors of ^44j derived the algebraic equation as the following limit of the Baker- 
Akhiezer system: 

F{P,Q)=0 s.t. F(^C,9str-g^)^{t;C) = ^?str ^ 0. (2.32) 

It was shown by the general arguments of topological recursions [10] that the algebraic 
equation has enough information to recover the perturbative corrections of the string 
theories. This means that we can reconstruct "perturbative Baker-Akhiezer functions" 
from the perturbative pieces of the algebraic equation: 

*^ilt(i;C) : F^^^(c,g,,,—^¥jUt;0 = g^tr ^ 0, (2.33) 
the leading behavior of which is given by one of the branches of the algebraic equation 

7^(^•)(t;C),F(c,7^(^■)(^;C)) =0: 

"^^pUt; C) ^ v^^\t; C) exp dCn^^\t; (')] +■■■, Qstr ^ 0, (2.34) 

asym L / 

where v^^\t;() is a proper fc-order vector- valued functionjf] Interestingly, these wave 
functions ^1/2-^ U = 2, ■ ■ ■ ,k) are completely decouple to each other in all-order per- 
turbation theory, since the algebraic equations are factorized into these irreducible pieces. 
Clearly, the exact Baker-Akhiezer function \l/(t;C) and the perturbative Baker-Akhiezer 
functions ^I/pert(^)C) different, but an asymptotic expansion of the exact Baker- 
Akhiezer function \l/(t;(^) in Qstr is given by a superposition of the perturbative Baker- 
Akhiezer functions \&p2rt(^i 0- 

^(t;C) ^ y^Cj^iJ^it-X), gstr^O, with (t,C)G^DcC2. (2.35) 

asym ' ' 
j 

This result suggests that the total wave function of the fractional superstring theory is a 
superposition of various perturbative minimal string theories [H] ! One of our motivations 



"'Each equation F'^^\P, Q) = corresponds to an eigenvalue of the Lax pair (P, Q). Then the upper 
index of the equations satisfies F^^-^^ (P, Q) = p(j) (P, Q) = F^j+^^ (P, Q) = 0. 

^ These perturbative Baker-Akhiezer functions also receive non-perturbative corrections by ZZ-branes 
of the corresponding perturbative string theory. The general form of these functions is universal and is 
expressed by theta functions on the spectral curve [42] (in other words, they are the r-functions of the 
corresponding integrable hierarchy). 
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in this paper is to further investigate this intriguing point from the non-perturbative point 
of views. 

In particular, from the perturbation theory, there is no way to fix the relative co- 
efficients {cj}j^i. This is the orthonormal-polynomial version of the non-perturbative 
ambiguity, and it is clear that these coefficients {cj}j^i cannot take arbitrary values, 
even though these parameters {cj}j^i are an analogy of the theta parameters in QCD. 
The correct values are obtained in the Z^-symmetric critical points [72| and we will show 
the fractional-superstring cases in Section 13.31 Interestingly, as is shown in [72] and as 
we will see in Section 13.31 determination of these coefficients {cj}j^i are equivalent to 
fixing the D-instanton fugacities. 



2.2.3 The multi-cut geometry and Stokes phenomena 

By an analytic continuation of (, the most dominant perturbative Baker- Akhiezer func- 
tions in the asymptotic expansion Eq. fl2.35p will change in a discontinuous way. Such 
discontinuities generally appear when one crosses the Stokes lines: 

SL,,, = {CGC;Re dC[n^^\t- (') ~ n^'\t- (')] = o}, (j, / = 1, 2, ■ ■ ■ , fc). (2.36) 

In particular, the Stokes lines can be interpreted as physical cuts of the eigenvalue distri- 
bution function p(A) of Eq. (IXSllFI if the exact Baker- Akhiezer function is given by the 
scaling orthonormal polynomials \E'orth(^; C) of Eq. (12. lip (and Eq. (I2.10p ). 

^{t; C) = ^orth(t; C) ^ v{t; C) exp \g;,l f dQ' niQ')] + • • • , 

asym L / 

g^tr -> with (t, C) (^^D C (2.37) 

Here 7^(C) is the scaling resolvent operator, NdxR{x) = g^^l d(^TZ(Q. From the viewpoint 
of the large N resolvent, Eq. (12. 8p . it is clear that the physical cuts (around ( — )■ oo) 
appear along the particular angles, C — oo x e*^"j3 

Xn = Xo + -^, Xo = -^il-e), (n = 0,l,2,--- ,fc-l). (2.38) 

However this consideration is not straightforward from the Baker- Akhiezer function sys- 
tem Eq. (I2.14P and implies non-trivial physical constraints on the coefficients {cj}j^i 
in Eq. (I2.35P and also on the Stokes phenomena which relate the coefficients in differ- 
ent regions D of the asymptotic expansions Eq. fl2.35p . This is the multi-cut boundary 
condition proposed in [72] . 

This consideration in the fractional-superstring critical points is quite interesting be- 
cause the perturbative minimal string theories described by the irreducible algebraic 
eqations, F^^\P, Q) = 0, (j = 1, 2, ■ ■ ■ , k) are completely decouple in the all-order per- 
turbation theory. Therefore, the discontinuous jumps of the wave function \E'orth(^; C) 
are the jumps between distinct perturbative string theories (i.e. not the jumps within 



^The physical cuts can be curved lines if the matrix-model potential includes complex coefficients. 
However the eigenvalue distribution function along the lines should be a real function. This requirement 
is guaranteed by the definition of the Stokes lines, Eq. (|2.36p (See [72]). 

^Note that there is a nontrivial rotation of angle by e, a; '-^ '^'^^^C: in Eq. (|2.12p . 
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the same perturbative string theory). Therefore, this means that different perturbative 
string theories are distributed in different regimes of ( and they are connected by Stokes 
phenomena! This idea was first introduced in [H] in order to add physical cuts which 
disappear in the spectral curves F{P, Q) = 0. A quantitative demonstration of this idea 
was then carried out in the symmetric critical points [72]. Here we now analyze what 
actually happens in the fractional-superstring critical points. 



2.3 Stokes phenomena in the multi-cut critical points 

The general frameworks of the Stokes phenomena in the multi-cut critical points are 
developed in Section 3 of [72] (a review of basics of Stokes phenomena with relevant 
references is in [87] or in Section 2 of [72]). Here we will not repeat the discussions but the 
relevant results are briefly summarized in order to fix the notations and conventions. Some 
new features appearing in the fractional-superstring critical points are also mentioned. 
For sake of simplicity, the discussion is restricted to p = 1. 

In the case of p = 1, the Lax operator P{t; () is generally given as 

P{t;d) = Ad + H{t), detA^O, (2.39) 

and therefore the differential equation system Eqs. f l2.14p is rewritten as the following 
ODE system of the Zakharov-Shabat eigenvalue problem [90|[9T]: 

<7str^^(t; C) = V{t; C) ^{t; ^ (c - H{t)^ ^{t; (); (2.40) 
9str-^^{t; = Qit- xl>(t; C) = Q{t; d) ^{t; C). (2.41) 

The Douglas (i.e. string) equation [T2] is now given by the following Zakharov-Shabat 
zero-curvature form of {V, Q): 

[P{t-d),Q{t-d)] =g,,Jk ^ [9strg-^-V{t-C),9strg^-Q{t;0] =0. (2.42) 

This system is referred to as A; x A; isomonodromy system, since this differential-equation 
system preserves the monodromy data and therefore Stokes data of Eq. (I2.4ip with respect 
to the flows of (t; (). 

The leading behavior of the ordinary differential equation, Eq. f l2.4ip . is given by 

9str^^^ = {A-X^-' + . . .)^it; C), 7 = r - 2. (2.43) 

Here r = g + 1 is the Poincare index of the essential singularity at ^ — ?■ oo. We consider 
two choices of the matrix A^ 

. _ J r ~ : w-'^/^-rotated models n — a- ^ ( i-i^ (o aa\ 

I r(,eai) ^ uj''^^ : real-potential models ' = ^'^^j=i{^ )■ (^-44) 



^In this paper, "~" indicates equality up to some similarity transformation. 
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As in [72], we also restrict ourselves to the cases of g.c.d. (A;, 7) = 1 in order to avoid 
complexity due to degeneracy of the exponents and appearance of the subdominant ex- 
ponents. Interestingly, this condition coincides with the following condition on dg^p in 

dg.p = g.c.d. ik,q-l) = l. (2.45) 

From the fractional-superstring viewpoint, this condition means the restriction to the 
critical points which are free from the x Zfc-orbifolding [73]. This kind of orbifold- 
ings was first introduced in the two-cut critical points of the two-matrix models [39] 
(i.e. Z2-orbifolding in the cases) in order to keep the correspondence with perturbative 
type minimal superstring theories of odd sequence {q,p G 2Z + 1). Therefore, the 
above coincidence suggests that this Z^ x Z^-orbif aiding procedure is naturally encoded 
in the non-perturbative Stokes-phenomenon structure and can be imposed by just com- 
pletely identifying the degenerate exponents. This is very interesting observation but the 
investigation along this direction is out of scope in this paper. 
The basic ingredients of Stokes phenomena are following: 



2.3.1 Asymptotic expansions 



This ODE system Eq. f l2.4ip has the k independent order-A; column vector solutions 
'^^■'\t; (), {j = 1,2, ■ ■ ■ , k), and they are often expressed as the following matrix: 



^{t;0^{^¥'\t-x),--- ,^^'\t-o). 



(2.46) 



The Stokes phenomena discussed in this paper originates from the asymptotic expansion 
around ( — i- 00. In term of the matrix notation Eq. f l2.46p . with a proper change of 
normalization of the solutions, the expansion can always be expressed as 



^asymit;C)^Yit;C)e±^''^'-'^^ 



n=l ^ 



X exp 



— ( V^o In C 

Qstr 



m,=— r,7^0 



(2.47) 



The leading exponents are written as 

^{t-X) = — 



C + 0{C 



C 00. 



(2.48) 



From the hermiticity of the matrix models, these matrix solutions are always real func- 
tions [12]. This basis of the solutions is called the matrix-model basis (or the T -basis). 
On the other hand, the Stokes phenomena are always discussed in the diagonal basis (or 
the Q-basis): 



asym 



{t; ^ Y{t; e^^^'-'^^ 

" Y^ty 



4 + E 



n=l 



X exp 



fi'str 



m=—r,^0 



— ^ asymi,^: C) ^) 



(2.49) 
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where the unitary matrix U is given as 

f/^, = 1 ^a-l+e/2)a-l+e/2)^ AU = U U^'Q) . (2.50) 

Therefore, the exponents are diagonal matrices: ^(t; C) = diag^^^ ((/^(■'^(t; C)), and the 
leading behavior is always the same as the exponents of the perturbative Baker- Akhiezer 
functions Eq. ( ICTll . 

V^'\t;0 - j' dC' n^'\t-X') = 0(C°), C ^ oo. (2.51) 
2.3.2 Stokes sectors 

The reconstruction of the exact solutions from these asymptotic expansions is achieved 
in specific angular domains called Stokes sectors Dn- This can be defined by Stokes lines 
Eq. (12. 36 p . especially their leading behavior of C — > oo: 



sl,v ^ {C G C; Re[(^^] - = O} 



C e e^'^U; 9[-^ = kn + ^{j + l-2 + e) ^^ ^ ^^.52) 



where is the leading exponent of ip''-'\t; () = (1/r) if^llc + - ■ ■ , and 7 = r — 2. Because 
of g.c.d. (/c,7) = 1, the minimal angular domains in between the Stokes lines, 6Dn, are 
given bjo 

^^n = I^(^^,^), (2.53) 

which are referred to as segments. Therefore, the Stokes sectors are defined as the angular 
domains which consist of {k + 1) different segments: 

^»-^(^.^). (2.54) 

which are also called fine Stokes sectors^ In each Stokes sector, one can uniquely recon- 
struct the exact (canonical) solutions \E'ri(t; () from the asymptotic expansion Eq. (I2.49p : 

$„(t;C) ^ 5asym(t;C), C ^ OO G (2.55) 

asym 

However, in general, these exact solutions have different normalizations: 

^n+i{t;0 = ^n{t;OSn. (2.56) 

This behavior of asymptotic expansion is called Stokes phenomena in this ODE system 
and the matrix Sn is called (fine) Stokes matrix. The total number of segments is 2rk 
and is equal to that of fine Stokes sectors. 



2.3.3 Profiles and Stokes multipliers 

Reading the positions of non-trivial components in the Stokes matrices, S'„ = {sn,i,j), is 
always a tedious problem, although there is a basic theorem which can be applied to 
every ODE system in principleEl] Here, however, we use more efficient way which was 

^In this paper, we use the fohowing notation of angular domain: D{a, b) = {a < arg(C) < b}. 
^°This choice of Stokes sectors is most fundamental. For other descendant Stokes sectors, see [72] , 
^^See Theorem 2 in [72], for example. 
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found in [72] with profile of dominant exponents. The profile of dominant exponents J' 
for {k,r;'j) is a sequence of integer numbers {j;,„} of the following type: 





<^2rfc-l 




J2rfc-l,l 


J2rfc-1,2 




j2rk-l,k 


J = 














Jl 






Jl,2 








Jo 




Jo,i 


J0,2 




jo,k 



(2.57) 



Each row vector J; corresponds to a segment of Eq. fl2.53p . and the numbers in 
Jl = [ ji,i I ji,2 I ■ ■ ■ I J«,A: ] describes the profile of dominant exponents in the segment 
5Di: 



Re < Re [^'I';'>C] < • ■ ■ < Re [^'I';''C] , C G 



(ii.fc) , 



(2.58) 



The profile is expressed as J'^^\^ in order to explicitly show the indices {k, r) and e = 0, 1. 



In this paper, 7 is always 7 = r — 2. As an example, a part of the profile J'l^^\^ (7-cut 
(P-iQ.) = (1)2) w^/^-rotated fractional-superstring critical point) is shown: 



(0) 



(7,3) 
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(5 


6) 
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(4 


6) 
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7) 


(2 


1) 


. (4 


5) 


(3 


6) 


(2 


7) 


1 



(2.59) 



Here we put the parentheses as in the profile to indicate the pairs which 
change the dominance in the next segment (i.e. there is a Stokes line of skj just above 
the segment). We also use the notation {i\j)i to show the segment J/ which the pair 
belongs to: € J/. Then the non-trivial claim (Theorem 4 in [72]) is 

si^ij is non-trivial component of Si 4^ (jK); ^ Ji- (2.60) 

The others off-diagonal components are zero, si^ij = 0. From the standard theorem, the 
diagonal components are always s/^j^j = 1 {i = 1,2, ■ ■ ■ ,k). 

2.3.4 Formula of profile components 

A simple formula for the integer numbers in the profile „ is also obtained in [72] (The- 
orem 3) for the w^/^-rotated models and the formula for the real-potential models is also 
obtained by a minor change 



(/ - 67) 



k — n + 1 



mi. 



mod k, 



(2.61) 



-•^•^Note that the first exponent ^^IC — ^ '^/'^C,^ becomes the largest when arg(C) = 1— 7- Also note 

kr 

that the effect of e is just a total shift of in Eq. (|2.52l) . 
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where mi is obtained by the Euchdean algorithm of krii + 7ml = 1. Interestingly, the 
difference of the formula for the Stokes phenomena between the cj^/^-rotated models and 
the real-potential models is just a shift of the angular coordinate and essentially they both 
have the same structure, even though these models are distinct from the viewpoint of 
perturbation theory. This implies that different perturbative non-critical string theories 
are simply unified in the strong-coupling non-critical M-theory description, as is in the 
critical 11-dimensional M theory. 



2.3.5 Three basic constraints 

There are basic three (or two) constraints on the Stokes multipliers. The original three 
constraints in the two-cut cases can be found in ^87j and also the extension to the multi- 
cut cases is in [72] with referring to the matrix-model results |13]. A major difference in 
the fractional-superstring critical points is that there is no symmetry (or there is only 
Z2 symmetry when k G 2Z). These constraints are also obtained in the similar way as 
in 



Z2-symmetry condition (only when k G 2Z): 

Sn+rk = T--2SnT-2, = 0, 1 , ■ ■ ■ , 2r - l) . (2.62) 



The hermiticity conditions: 

S: = Ar(i-) 5(2Li).-n r-(^-^) A, = 0, 1, ■ ■ ■ , 2r A; - 1) , (2.63) 

with the matrix A, Ajj = 6i^k-j+i- Note that the unitary matrix U of Eq. fl2.50p satisfies 
U"^ = (— l)'^Ar^^^'^''. In terms of the elements, they are expressed as 

= -S-k-n,2-e~i,2-e-j, (jK)n ^ J„. (2-64) 



The monodromy free condition: 

Sq ■ Si - ■ ■ S2rk-1 = Ik- 
la the case of A; G 2Z, this is also expressed as 

(>S'o ■ ■ ■ S'r.fc_l)r 2 = 

In the case of G 2Z + 1, this is also expressed as 

s* = Ar(^-^)5Ar(^-^), 

with S = S_k+i_^^ S_ 



(2.65) 



(2.66) 



(2.67) 
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See also Section 2 in (75] since the convention used in the matrix models is different. 
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3 The multi-cut boundary conditions (BC) 



In this section, we evaluate the multi-cut boundary conditions in the fractional-superstring 
critical points. The multi-cut boundary condition is first proposed in [72] (Definition 10) 
and also mentioned around Eq. fl2.38p in this paper. 

One of the main differences from the symmetric cases |72j is that there are a huge 
number of Stokes multipliers^ It is because the smallest Poincare index r is r = 3 
(there is no (p, q) = (1, 1) critical point) and there is no symmetry. This causes a 
predicament in direct analysis of the algebraic equations of Stokes multipliers. Therefore, 
we should invent some systematic way to analyze these large systems. 

An important property of the formula for the profile components Eq. f l2.6ip is the 
following vertical shift in the profile: 

jSa,n=4n+«^l- (3-1) 

This implies the following fact: 

si^ij is non-trivial (i.e. G Ji) 

Sl+2a,i+ami,j+ami IS UOU-trivial (i.e. "'(j + ami\i + ami)i+2a e Jl+2a)- (3.2) 

Therefore, once we know the existence of non-trivial Stokes multipliers of si^i^j {I = 
0,1,2,-- - ,2r— 1), then one can guarantee the existence of the other Stokes multipliers 
by using the above shift rule. This rule is referred to as 2a shift rule to indicate the range 
2a of making the shift. 

In the following discussion, it will become clear that the fundamental set of Stokes 

f ^ 2r— 1 

matrices in the above sense is \Snf^_^ , therefore essentially a single symmetric Stokes 
matriceJ^ Sq^"^\ Here the symmetric Stokes matrices are defined as 

^2rn ^ = S2rnS2rn+l " " " S2r{n+l)-l, ('^ = 0, 1, 2, ■ ■ ■ , A; — l) . (3.3) 

r ^ k—1 

Therefore, we will start with the symmetric Stokes sectors \D2rn} in the following 
discussion of the multi-cut boundary conditions. 



3.1 The multi-cut BC recursion equations 
3.1.1 The boundary conditions and the shift rule 

The scaling orthonormal polynomial ^orth(^; C) of Eq. (12.111) is generally expressed by the 
canonical solution \l/2rn(^; C) of the symmetric Stokes sectors -02™ as 

/^(2r'n)\ 

^orth(t; C) = ^2rn{t; C) ^^"""^ = 

^^For example, there are seven times more than the multipUers shown in the profile, Eq. (j2.59p . in the 
7-cut r = 3 case. 

^^The name comes fi'om the symmetric Stokes matrices in the Z^^-symmetric critical points |72) . al- 
though there is no Zj,-symmetry in the fractional-superstring critical points. 



(2rn) 



(2rn) I 



(3.4) 
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for each index n = 0,1,2, - ■ ■ ,k — 1. Therefore, the vectors X*^^*"") are related by sym- 
metric Stokes matrices of Eq. fl3.3p as 



x^(2rn) _ c(sym) y(2r(n+l)) 



(n = 0, 1,2, ■■■ 



(3.5) 



Then the multi-cut boundary condition is imposed on the vectors |X*^^''"')}^_^. By fol- 
lowing the procedure discussed in [72], we read the boundary conditions as follows: Here 
we first consider an example of 7-cut r = 3 cu^^^-rotated case. The profile in Dq and the 
corresponding vector X^^^ are given as 
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xy ^ 



v 



X 



(0) 

5 

(0) 
6 

(0) 



(3.6) 



Note that "2^" means xf^ 



and "4 A" means a; 4"'' 7^ 0. In the middle of the above 
profile, there is a horizontal line where "1^" "4^" change their dominance, which is 
interpreted as a physical cut of the resolvent operator. In principle we should perform 



r ^ k—1 

the same procedure for all the symmetric Stokes sectors |-D2rn}„_Q to obtain their vec- 



fc-i 



are 



tors however the boundary conditions for all the vectors {^^^'^"^}^_q 

automatically obtained from the first vector X^^^ in the first symmetric Stokes sector Dq 
because the multi-cut boundary conditions for the vectors are related to each 

other by the 2r shift rule of Eq. liS. 



I ^ / (2rn) 

7^ I or x!;- 







(2rn+2r) 
j+rmi 



7^0 



or 



(2rn+2r) ^ ^ 



(3.7) 



For example, one can directly check the other vectors {-^'•^"■' }^_-,^ of this case {k = 7, 
3, nil = 1) are given as 
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7^0 



X 



X 



X. 
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(12) 



(18) 



(18) 



X 



(18) 



7^0 








„(18) 



7^0 



V 4"^ / 



(3.^ 



and they are all consistent with Eq. (13. 7p . This is because we read the multi-cut boundary 
conditions from the profiles which are again obtained by using the 2r shift rule. Therefore, 
only the multi-cut boundary condition equation, Eq. (13.51) . of n = is essential and then 
the other equations are generated by using the 2r shift rule, Eq. (13.11) . This procedure 
reduces a great deal of the complexity. Note that this consideration does not mean there 
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is the Zfc symmetry in the critical points. An analogy of this consideration is "a symmetry 
of equation of motioij^." which does not necessarily mean a symmetry of the solutions. 

It is not difficult to read the boundary condition for general k with r = 3, which will 
be used in the practical analysis. We put k = 2rm + l (1 < / < 2r) and there are 2r (= 6) 



different cases 




= 1,2,- 




, 6, respectively): 










1. ^ (1^ 


X, 


X 4^ 


X, 


x,7^,-- 


■ , X, X, 


(3m + 1)^, 




); 




2. o (1^ 


X, 


x,4^ 


X, 


x,7^,-- 


■ , X, X, 


(3m + 1)^, 


X, 


); 




3. ^ (1^ 


X, 


x,4^ 


X, 


x,7^■■ 


■ , X, X, 


(3m + 1)^, 


X, 


X, ); 




4. o (1^ 


X, 


X 4^ 


X, 


x,7^,-- 


■ , X, X, 


(3m + 1)^, 


X, 


X, (3m + 4)^, 


-,-); 


5. ^ (1^ 


X, 


x,4^, 


X, 


x,7^,-- 


■ , X , X , 


(3m + 1)^, 


X, 


X, {3m + 


-,--); 


6. ^ (1^ 


X, 


x,4^ 


X, 


x,7^,-- 


■ , X , X , 


(3m + 1)^, 


X, 


X, (3m + 4)^, 


' 5 



(3.9) 



To save the space, we used a simpler notation: "1^" means Xi ^ 0, "x" on the i- 
th component means xf^ = 0, and "— " means that there is no constraint on these 

components. For example, the vector X^'^^ of Eq. (13. 6p is expressed as (1^, x, x, 4^, 

— ). Here we introduced a notation of the indices with tildes, (m, /). It is because this 2r 
classification is reduced to be a simpler r classification, k = rm + / (1 < / < r), in the 
final unified form of the multi-cut BC recursion equations!^ 

Finally we mention the meaning of the general 2a shift rule for the boundary con- 
dition. We used the 2r shift rule because the multi-cut boundary condition Eq. (13. 6p 
(generally mentioned around Eq. (12.381) ) only respects the 2r shift rule. Therefore, one 
may wonder if we use the other 2a shift rules on these multi-cut boundary conditions. 
The answer to this question is following: The 2a shift rules (a = 1,2, ■■■ ,r — 1) for the 
multi-cut boundary conditions generate the complementary boundary conditions. The 
complementary boundary condition on the Stokes phenomena was first discussed in the 
Zfc symmetric cases ^72j. The physical meaning of the complementary boundary condi- 
tions is still unclear but these constraints provide strong equations which help us obtain 
explicit expressions of the Stokes multipliers even in the general {k, r) isomonodromy 
systems!^ In fact, one can check that the 2a shift rule (a = 1, 2) for the the multi- 
cut boundary condition Eq. (13. 6p results in the following two kinds of complementary 

^^Maybe of non-critical M theory. 

^''Therefore, 1. O 4., 2. 5., 3. O 6. in Eq. (13.91) have the same form of the multi-cut BC recursions. 

-"^^One may expect that it is a constraint on operators in non-critical M theory which is analogous to 
"degenerate operators in conformal field theory". Also note that, if the indices (fc,r) are small integers, 
then the complementary boundary conditions are exactly satisfied. See |72j for explicit analysis of some 
concrete examples. 



18 



boundary conditions: 



2 shift: Do D 



4 shift: L>4 3 
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(3.10) 



(3.11: 



Therefore, in general, each of these 2a complementary boundary conditions (a = 1, 2, • • • , r- 
1) requires that there is a solution of the Baker- Akhiezer function system Eq. f l2.14p which 
has the physical cuts along the k directions, C, ^ 00 x e^^^ , of 



Xn = XO + 



2im 



Xo 



+ 



27ra 



(n = 0,l,2,--- ,fc-l). 



(3.12) 



k k rk 

In the following, it is also convenient not to distinguish the multi-cut boundary condition 
and the 2a complementary boundary conditions (defined by Eq. fl3.12p ). therefore, they 
are equally called the multi-cut 2a boundary condition. In the a = case, we skip 2a in 
the name. 



3.1.2 The theta parametrization 

In the analysis of the Z^-symmetric critical points [72j , it was observed that the multi-cut 
boundary condition equation Eq. (13. 5 p with the multi-cut boundary condition Eq. (I2.38P 
only includes a particular sequence of Stokes multipliers. Motivated by this observation 
given in [72], here we also introduce the following parametrization of the fine Stokes 
multipliers: 

^r+^^^^S .(.) .(e) 

2r+2a-n+l+er, J2r+2a-n+l+er,k-n+l ' 32r+2a-n+l+er,k-n 

(n = l,2,-- - ,A;-l;a = 0, 1,-- - ,rA;-l). (3.13) 

It is much easier to catch the meaning of the formula by drawing it on the profile. Figure 
[2] is an example of the profile for the 11-cut (r = 3) w^^^-rotated models. 

In fact, the multi-cut BC recursion equations Eq. (13.51) with the multi-cut boundary 
condition only include a subset of the Stokes multipliers, {^n^'^^'' }i<„<;fc„i! we expected. 
Importantly, the 2a shift rule acts on the parameter ^i^*^^^ as 

2a shift: ^ 0i2r6+2a)^ (3_^4) 



19 



9 


8 


10 


7 


11 


6 


1 


5 


2 


4 


3 


8 


9 


7 


10 


6 


11 


5 


1 


4 


2 


3 


8 


7 


9 


6 


10 


5 


11 


4 


1 


3 


2 


7 


8 


6 


9 


5 


10 


4 


11 


3 


1 


2 


7 


6 


8 


5 


9 


4 


10 


3 


11 


2 


1 


6 


7 


5 


8 


4 


9 


3 


10 


2 


11 


1 



(a) 











93*'^ 


9/°' 


A (8) 
f 10 








9/) 










93^^) 


93*"^ 


0^(64) 


A (6) 
f 10 


A (4) 




94^°^ 














93(^4' 




62) 


A (4) 

fio 








64) 


Q2 


62) 





(b) 



Figure 2: The theta parametrization of fine Stokes multipliers in the 11-cut (r = 3) w^/^-rotated 
models. Here only the fundamental domain IJj^q J; of the profile is shown, a) The profile of dominant 
exponents with classification by different 3 {— r) colors. Here we skipped the parenthesis notation 
but each pair of two same-colored boxes next to each other is the parenthesis pair b) The theta 

parametrization 0n'^^ is shown by using the correspondence s;,ij -f-)- S Ji- The colors are drawn in 
accordance with a modulo r. One can see the tilting pattern of the parametrization. 

therefore the multi-cut BC recursion equations Eq. (13. 5 p with the multi-cut 2a boundary 
condition (a = 1, 2, ■ ■ ■ , r — 1) include a complementary set of the Stokes multipliers, 

r^(2rfe+2a)-i 0<fe<fc-l 

l*^" /i<„<fc_r 

Before we discuss the details of the multi-cut BC recursion equations, basic properties 
of the theta parameters are summarized: 

• The Z2-symmetry condition {k G 2Z) Eq. f l2.62p is expressed as 

Qi2a) ^ Qi2a+rk) = 1, 2, ■ ■ ■ , A; - 1; a = 0, 1, 2, ■ ■ ■ , rfc - l) . (3.15) 

• The hermiticity condition Eq. (I2.63P is expressed as 



(n = 1,2,--- a = 0,l,2,--- ,rfc-l). (3.16) 



For later convenience, we extend the parameter as 

^(''^^ = -1, ^f) = +l, (a = 0,l,2,---,rA;-l), (3.17) 
which is consistent with the hermiticity constraint. 



3.1.3 The multi-cut BC recursion equations 

In this section, we discuss the multi-cut boundary condition equation Eq. (13. 5 p with the 
multi-cut 2a boundary condition Eq. (I3.12p . The derivation is given in Appendix |Al and 
therefore we here only show the result: 

Theorem 1 (The multi-cut BC recursion equations) In the fractional- superstring 
critical points of an arbitrary {k,r) with g.c.d. (/c,7) = 1 ('-/ = r — 2), the multi-cut bound- 
ary condition equations Eq. I{3.5\) with the multi-cut 2a boundary condition Eq. h3.12^) of 
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a = 0, 1,-- - ,r — 1 result in the following recursion equations: 



J", 



i2s) 



of 



0<r(n-l)+6<fc,neZ 

r/ + a + e^; (O < / < A; - 1, l<fe<r), 



where the sequence of numbers {nn^^^^^ ore all non-zero and satisfy 



(3.18) 



,,(2.) _ (2s+2r6) 
i/n ~ yn-b ' 



y{2s) _ (2S) 

iln yn+k^ 



[n, a,b E Zj. 



(3.19) 



Therefore, there are r distinct sequences {l/n^"''}^_-^ corresponding to each boundary 
dition of 2a. 



con- 



(2s) 



of 6 



For each boundary condition (of 2a) there are r recursive equations J^i 
1, 2, • ■ ■ , r. It is suggestive to see the positions of the Stokes multiphers (included in each 
equation J^^^^*-* = 0) in the profile. Here we do not show exphcitly, but one can see that 
the trajectory is along the line vertical to the trajectories of colors in Figure [2j 

Although the formula Eq. 03.181) is simply written in a unified way, the concrete form 
of each equation system depends on k of modulo r: k = rm + /. Here we show each of 
them. First of all, we introduce a new parameter 



a, 



(2a) 



^{2r{n-l)+2a) 
r(n-l)+f) ' 



m—n+l,b 

and then the multi-cut BC recursion equations are 



k-1 



(3.20) 



J", 



(23) 



^(2«) i2s) _ Q 



(3.21) 



The concrete form of each equation is shown in Table [TJ Here we also show the hermiticity 
condition Eq. fl3.16p in terms of the sigma parameter: 



cr, 



(2s) 



m— n+1,. 



a 



(-2r(m+2)-2s) 
n-l,{l-b) 
(-2r(m+l)-2s) 
n,{r+l-b) 



: 1 < 6 < /, (ra = 1,2, ■ ■ ■ ,m + l) 
: I < b < r, (ri = 1, 2, ■ ■ ■ , m) 



(3.22) 



Therefore, the hermiticity condition relates equations of the five cases (I, II, HI, IV, V) 
as: 



lol, H^IV, IH^IH, VoV. 
These are consistent with Eq. f l3.2ip (also see Tabled]). 



(3.23) 



3.2 Relation to discrete Hirota dynamics 

Here we point out an intriguing connection to the structure of quantum integrable models. 
A non-trivial connection between Stokes phenomena in ordinary differential equations 
and T-systems of quantum integrable systems is found in a special kind of Schrodinger 
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cases 


1 < 6 < r 




length 


I. 


l<b<l 


7Ti-|-l 

^(2^) y(2s) _ g 
^'^■m-n+l,byn — 

n=l 


m + 1 


II. 


h = l{^r) 


771 

+ V a^^"^ - 

i/m+l ^ /_^"m-n+l,b iln ~ ^ 
n=l 


m + 1 


III. 


I <b < r 


m 
n=l 


m 


IV. 


h = r{^l) 


m 

yo ^ 2-^ '-'m-n+l,b y-n — ^ 
n=l 


m + 1 


V. 


b = r{=l) 


m 

(2s) V- (2s) (2s) _ 
n=l 


m + 2 



Table 1: The concrete forms of Eq. p.lSp (and then Eq. p.2ip ) is shown with k = rm + 1, {1 < I < r). 
The case of V is the degenerate case of II and IV (i.e. r = I). Possible cases are (II, III, IV) for I = 1, (I, 
II, III, IV) for 1 < I < r, and (I, V) for I = r. The "length" means the number of terms in the equation 



equations [921193] . which is called the ODE/IM correspondence^ In this section, we 
provide a new generalization of this correspondence to general ODE systems appearing 
in the isomonodromy theory. 

In general, the Stokes multipliers are distributed in the Stokes matrices 
in a very complicated way. However, if one looks at the general multi-cut BC equa- 
tions Eq. fl3.18p . one notices that the Stokes multipliers are re-organized in the following 
suggestive form: For example, the equation of IV in Table [1] is expressed as 



y(2s-20 ^ ^(2s) ^(2s) ^ (s = 0, 1, 2, ■ ■ ■ , rA; - 1) , 



(3.24) 



with 



y(2^) = 



/ (2s) (2s) 

/ 1 / O 



$(2^) 



1 





1 



(2s) 
m—1 






I'm \ 






V 



(3.25) 



/ 



Here {''n^'^'' }^gg are related to the fine Stokes multipliers {(^nb^nz 



(2s)- 



-(2s) 



(2s) 

m— ri+l,r ' 



( 



n 



1,2, 



m 



(3.26) 



This is the same form of the recursive equations found in the study of Stokes phenom- 
ena for the Schrodinger-equation systems of the ODE/IM correspondence (See [93J for 
example). Therefore, according to the discussion in [93], each sequence 



""^^The authors would like to thank Toshio Nakatsu and Ryo Suzuki for drawing our attention to these 
interesting papers and the correspondence. 
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(6 = 1, 2, ■ ■ ■ , r) satisfies discrete Hirota equations [HI] (or a T-system) of quantum inte- 
grable systemso An important difference from the original correspondence is that in the 
k X k isomonodromy ODE system of the Poincare index r, there are r different recursive 
systems of Stokes multipliers which result in r different T-systems. That is, the system 
is described by multiple T-systems of quantum integrable systems. 

We note that, in our calculation, there is no explicit appearance of spectral parameters 
of quantum integrable systems. However the spectral parameters appear by assuming 
analyticity of the index 2s. That is, we rewrite the index s as s = ru + a {0 < a < r — 1) 
and express them as follows: 

r(2^) = t4u; a), ^ = yo{u + n; a), Y^'^^ = Y{u; a). (3.27) 

Then the claim is that analytic continuation of the index u realizes the spectral parameters 
of the T-systems. This is understood as continuous deformation of the solutions for Stokes 
multipliers. Because of the isomonodromy property, therefore, this spectral parameter u 
commutes with the original spectral parameter ( of the resolvent operator: 



d_ dj 

du dC 



^{t-X;u) = 0. (3.28) 



In this sense, this is a new spectral parameter which describes the strong-coupling dy- 
namics of non-critical string theory, possibly, non-critical M theory. 

3.2.1 Hirota dynamics for non-critical string/M theory 

Here we skip "a" which does not change in each integrable model. The direct consequence 
of [93j is following: T-system appears as the algebraic relations encoded in the recursion 
equation Eq. fl3.24p and Eq. fl3.25p . i.e. in the algebraic relations among components 

(u) = {(j)n,i,j{u)) ^^: 

$„(m) = ^{u) + l) ■ ■■^{u + n- 2) ^{u + n- l). (3.29) 

The discrete Hirota equation (i.e. T-system) which is relevant to the algebraic structure 
is that of the integrable models with the quantum symmetry [/^(A^Li) of g'^ = 1 [93] : 

Ta,s{u + 1) Ta,s{u - 1) = Ta,s+l{u) Ta,s-l{u) + Ta+l,s{u) Ta-l,s{u), (3.30) 

with the following boundary condition: 

T.nAu) = Tm+nAu) =0 (n G N, s G Z) , (3.31) 
Ta^s+k-m{u) = Q (O < a < m, < s < m), (3.32) 
Tafl{u) = Tq^u) = I (0<a<m, sgZ), Ta^u) = Ta^s+k{u) , (3.33) 

which is also shown in Figure El Note that T-functions also satisfy periodicity in the 
spectral parameter u: 

Ta,b{u + 2k) =Ta,b{u). (3.34) 
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For references of Hirota equations, see [SS] for example. 
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The correspondence of T-functions with the components of the is given in the 

following way [93j0 

Tn{u) = (-l)"+iT„„i(2u + n), (t)n,iM = TiA'^u + n). (3.35) 

Notes on the boundary conditions are in order: 

1. Generally it is known that general T-functions of the above T-system can be ex- 
pressed in terms of the vertical/horizontal inputs by the Bazhanov-Reshetikhin 
formula [96] 

Ta,s{u) = det Ta+i-j,i{u + s + 1 - i - j), (3.36) 

l<i,j<s 

= det Ti^s+i-jiu + a + 1 - i - j), (3.37) 

l<i j<a 

with the boundary condition Eq. (13.331) . 

2. One can directly solve Eq. (I3.29P and check that the coefficients {(pn,i,iiu)}'^_^ 
are expressed by {r„(u)}™^^ and the expression is equivalent to the Bazhanov- 
Reshetikhin formula of Eq. (13.361) (and the identification Eq. (I3.35p ) with the fol- 
lowing extension of an index: 

ro(u) = -l, r_„(n) = w„,(n) = {neN). (3.38) 

In particular, Eq. (13.361) results in Eq. (I3.3ip and Eq. (I3.33p . 

3. The cyclic condition ^k{u) = Im (which results from the monondromy free con- 
dition Eq. (12.651) in our case) gives the following constraints on the components 

{<^n,l,l(M)}„g^: 

4>n,l,l = (t)n+kXl (neZ), (j)k~m+n,l,l = Sn,m (l < U < m) . (3.39) 

Therefore, by using Eq. (13.370 . one can show Eq. (13.320 and Eq. (I3.33p . 



3.3 Explicit solutions of Stokes multipliers 

Here we show several explicit solutions for the Stokes multipliers which generalize the 
discrete solutions found in the Z^-symmetric critical points [72]. Solutions shown here 
are therefore special solutions and further analysis for the complete solutions (with using 
the above T-systems and their Bathe ansatz equations) should be studied in another 
publication. In particular, we impose the following "Z^-symmetry condition" : 

Sn+2r = T-'SrX (^ = 0, 1, ■ ■ ■ , 2rfc - l) 

^ 0^2s) _ ^^2s+2r) (l<n< s = 0,2,--- ,rA:-l). (3.40) 

Note that this constraint on the Stokes multipliers does not necessarily mean existence 
of the Zfc symmetry of the system, although the monodromy free condition Eq. (I2.65P 
becomes 

[St^'^r-'f = Ik. (3.41) 
^-'^The expression for more general components is also shown in [53] (i.e. Theorem 3). 
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i a 

; m 



k-m 



Figure 3: The boundary condition for the T-system with = 1- The lattice points 

exphcitly shown in the figure are the support of the T- function, To,s('w)- In particular, To_s(w) = Tafl(u) = 
1 (s G Z,0 < a < to), and also there is a periodicity Ta^siu) — Ta^s+k{u)- The example shown here is 
fc = 14 and m = 4 (i.e. r = 3). 



From the viewpoint of the discrete Hirota equation, the Zfc-symmetric constraint is equiv- 
alent to dropping the spectral-parameter dependence and the Hirota equation then be- 
comes the system for characters of a corresponding Lie group. Consequently, the Stokes 
multipliers {t„(m)}^__|^ of Eq. (13. 25^ are given by the characters of the anti-symmetry 
representations, Xn,i{g)- 

t4u) = {-ir+'T^,,{u + 2n) = {-l)-+\^^,{g), (3.42) 

with a proper group element g G GL(m). The constraints on the group element g is 
from the algebraic relation of the T-system boundary conditions, however, we know the 
solutions since the system is the same as this Z^-symmetric situation [72]. That is, the 
group elements are given by distinct roots of unity, jcj'^J characters are given 

by symmetric polynomials: 

XuAg) = Sym[n, {co-^jZ,] , 9 ^ diag(a;"S ■ ■ ■ (3.43) 



This is the connection to the discrete solutions found in 

Clearly this is not the end of story because there are now r different T-systems, 

Tn{u) rn{u;a) (a = 0, 1, 2, • • • , r - l) , (3.44) 

and they are related by the original monodromy free condition Eq. (12.651) . This condition 
can be solved by showing explicit k distinct eigenvectors of the symmetric Stokes matrix 
^(&ym) jY2j_ fjj-^jg condition results in the constraints on the charges (or group elements) 
g{a) which are classified by configurations of the avalanches \72\ . 



3.3.1 The case of k = rm + 1 

In this case, we have three kinds of equations (II, III and IV in Table [1]). Therefore, we 
have the following recursion-equation systems 

II: Y{u + l;a) = ^{u;a)Y{u;a), (a = 0, 1, ■ ■ ■ , r - l), 
III: = Vb{u;a)Y{u;a), (l < 6 < r; a = 0, 1, ■ ■ ■ , r — l), 
IV: Y{u-l;a) = <!?{u;a)Y{u;a), (a = 0, 1, ■ ■ ■ , r - l) . (3.45) 
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The matrix a) is defined by Eq. fl3.25p and Eq. (13.261) . and the others are defined as 



a) 

Vb{u; a) 
We here requirj^ 






1 







n 
u 


n 
u 


1 

i 


























{2s) 
-^m,l 


(2.) 
m— 1,1 


^(2.) 
'^m-2,1 




X2s) (2s) 
rn,b 1 ^m-l,fe' 


(2s) 
'^m-2,6' ■ 


' ^2,6 ' 


^(2s) 








\ 





1 










^(2s) 





\7) {l<h<r). 



$(u;a) = <l> ^(M+l;a), a) = 0, 



(3.46) 



(3.47) 



then the recursion equations is given by a single equation Eq. (13.241) discussed in Section 
13. 2[ which means that there are r Hirota equations. In terms of the components, these 
Stokes multipliers are given a@ 



II : 

III : 

IV : 



cr, 



cr, 



cr. 



{2s) 
m,l 

(2s) _ n 

{2s) 

m—n+l,r 



{2s) 
m— 71,1 



r, 



(2.) 



(2s) ^ "m-n,L J2s) 

' m I m 



n 



1,2, 



, m 



-1). 



(l<6<r;?2 = l,2,---, m) , 
f^), (n=l,2,---,m). 



(3.49) 



As special solutions to the Hirota equation for {r„(M; a)}°<^<^ ^ (discussed in the begin- 
ning of this section), we assign the following discrete solutions 



Tn{u- a) = (-l)"+^Sym n, {u;"^"^ j^^^ (a = 0, 1, ■ ■ ■ , r - l) . 



(3.50) 



Once we derive this formula, the procedure to obtain solutions is the same as the Z^- 
symmetric cases [72]. Here we explicitly checked the constraints by Mathematica, and 
therefore, we make the following conjecture for the conditions on the exponents of u in 



^^Note that there are other solutions which does not satisfy this condition. For example, we can also 
consider recursion equation of type III and define T-system of Uq{Al^^_2). Here however, we only focus 
on special solutions which are relatively simpler. 

^•^Note that the inverse of the matrix is given as 



/ 






1 

A2s) 








\ 







1 



(3.48) 
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Eq. em-. {nrYSt' 



(2a) I (2a) 



mod k 7^ j; < a < r — l) 

k\l + {-l)^ 



(3.51) 



mod k 



l<j<m 
0<b<r-l 



V-J /l<7<m 1 



l<jr<m 

Some notes are in order: 



n 



{2{r-a)) ■ 



l<j<m 



(1 < z < m; < a < r - 1). 
(0 < a < r - 1). 



(3.52) 
(3.53) 



In order to satisfy the monodromy free condition Eq. (I3.4ip . one should show that 
the symmetric Stokes matrix Sq^^™^ is diagonahzable. With the exphcit expres- 
sion of the discrete solution Eq. f l3.50p . one can explicitly construct rm {= k — 1) 
eigenvectors X{a;n^^°'^) 



^(sym)p_l ^2a]^ ^ ^nf^^ ^2a]^ (l < J < m; < a < r - l) . (3.54) 

A brief summary of the procedure is following: Each T-system of |T„(m; a)}™__^ 
(for a fixed a = 0, 1, ■ ■ ■ , r — 1) admits m different vectors Y{u; a): 



Y{u; a) = y{u; a; nf""^) = 



/^-(«+l)nf'\ 
-(«+2)nf''' 



(l < j < m) , 



(3.55) 



which satisfy 



$(«; a) y{u; a; nf^) = y[u- a; n'-'""). 



{2a) 



(2a), 



(3.56) 



By using the recursion equation Eq. (13. 5p . one can find the explicit expression of 
each eigenvector X{a] n^^""^) in terms of y{u] a; n^^""^): 



.,.if''^) = M{u-a)y{u-a-nf^), 



with a k X m matrix M{u; a) which satisfies 



(3.57) 



(3.58) 



{2a)x -I 0<a<r--l 



then one obtains the explicit expression for rm vectors {X(a; Uj )}i<j<„ 



. Since 

each multi-cut 2a boundary condition is distinct to each other, the vectors are also 
distinct if and only if the condition Eq. fl3.5ip is satisfied. 



Although the rm {= k — 1) distinct eigenvectors are explicitly constructed in the 
above procedure, there remains one missing eigenvector. Noting that det [5'o'^'^™''r^^] 
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(—1)'^ ^, one can conclude that the eigenvalue associated with this vector is given 
by w"" of 



no-- E + -odfc, (3.59) 

Cfe<r-1 

and the eigenvector is given as 

-^K) = J]c,M(M;a)3^(u;a;no). (3.60) 



l<j<m 
0<fe<r-l 



r-1 



a=0 



It is clear and the constraint Eq. ( 13.52^ is sufficient. The necessity is non-trivial 
but we explicitly check it by Mathematica for several cases. 

• These Stokes multipliers {^niu; a)^-^^^^^ should be subjected to the hermiticity 
condition Eq. fl3.22p . This constraint results in the constraint Eq. f l3.53p . 

3.3.2 The case of = rm + /, (1 < / < r) 

In this case, we have four kinds of recursive equations (I, II, III and IV in Tabled]), we 
also perform the same procedure shown above. Here we only show the identification: 



_(2s) ^(2s) (2s) 
^ mfi -L ^ m—n,b 'Tn 

~W) ~ i'is) ' (2s) ~ (2s) 



(ra = 1,2,-- ■ ,m- 1; 1 < 6 < /), 



II ■ a^^'^ - ^ a^^'^ - ^ (n-l2 ■■■ m - ll 

• "m,l — (2s)' ^m-n,l - (2s) l,"- - i, Z, ,171 I), 

Tni Tm 

III: aiy=0 {1< b<r;n = 1,2,- ■■ ,m), 

IV: ai^t+i,r = rt\ (n = 1, 2, ■ ■ ■ , m) , (3.61) 

then we have a single recursion equation Eq. fl3.24p which again results in r different T- 
systems. Therefore they are given by the discrete solution expression Eq. f l3.50p . However, 
the conditions for the indices are not straightforward in this case, and therefore, we show 
special solutions of simple cases. 

The case of {k,r) = (5,3) (i.e. (m, /) = (1,2)) 
We found two types of solutions. First we put 

: free, a^^f = (« = 1, 2), n = nf\ n = nf ^ = n^^K (3.62) 

Then the constraints are given as 

1. n^n mod /c, r^i = n + n, Ti2 = —n — 2n mod k, 

2. n = n mod k, rji + ri2 = —n. (3.63) 



The case of (A;, r) = (14,3) (i.e. (m, /) = (4,2)) 
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We found the following solutions: First we put 



cr, 



(2a) 
n,l 



(n = 0,1,2,--- ,m; < a < r- 1), 



nj = nf''^ (l < J < m; < a < r - l). (3.64) 



"J — '"J 

Then the constraints are given as 

m ^ 

'^nj = -, rij^k. (3.65) 
3.3.3 The case of A; = r(m + 1) 

In this case, we have two kinds of recursive equations (I and V in Table [1]). we also 
perform the same procedure shown above. Here we only show the identification: 



I: cr^y = (1 <6<r; n = 1,2,--- ,m + l) 



V: cjtt+i,r = r!^'\ (n = l,2,---,m), ri':\ = 1, (3.66) 

then we have a single recursion equation Eq. fl3.24p which again results in r different 
T-systems (i.e. each T-system is Uq{A^m) of = 1). To avoid confusion, we here show 
the boundary conditions which are slightly different from the other cases: 

T_„,,(u) = Tm+i+nAu) = {neN,seZ), (3.67) 
Ta,s+k-m^i{u) = (0<a<m + l, 0<s<m + l), (3.68) 
Ta,o{u) = To,s{u) = Tra+iM = 1 (0 < a < m + 1 , s G Z) , (3.69) 

Therefore they are given by the discrete solution expression Eq. fl3.50p . 

r„(«; a) = (-l)"+iSym [n, {uj^^"^^' (a = 0, 1, ■ ■ ■ , r - l) , (3.70) 
with the following conditions on the indices: 

1. nf"^ ^ nf"^ mod k (^ 7^ j; < a < r - l) . (3.71) 

m+l , I -1 \ 

2. = ^^'"^+^ "^ ^ ^ mod A; (0<a<r-l). (3.72) 

3. {"fV3^„={»f"''}.s,s„ (0<a<r-l). (3^73) 
Note that the fractional number in Eq. f l3.72p is always an integer. 

4 Conclusion and discussions 

In this paper, we analyzed Stokes phenomena in the k x k isomonodromy systems with 
an arbitrary Poincare index r, especially which correspond to the fractional-superstring 
(or parafermionic-string) multi-critical points {p,q) = (l,r — 1) in the multi-cut two- 
matrix models. Throughout the analysis, the multi-cut boundary conditions proposed 
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in [72] turn out to be a key concept in order to uncover the underling Hirota equa- 
tions (T-systems) and also a powerful tool in order to show explicit expressions of the 
Stokes multipliers in quite general class of {k, r). Although we here focused on fractional- 
superstring critical points of the multi-cut two- matrix models (i.e. 7 = r — 2), according 
to the procedure of imposing the multi-cut boundary condition discussed in Appendix 
Rl it is clear that the multi-cut BC recursion equations in Theorem [1] keep the universal 
form in the general critical points (i.e. arbitrary 7) of the multi-cut two-matrix models. 

It should be noted that the basic observables in this paper, Stokes multipliers, are 
integration constants of string equations (i.e. non-perturbative ambiguities) and are in- 
formation completely out of the perturbative framework of string theory. Therefore, it 
is quite striking to show existence of a rigid mathematical structure like quantum inte- 
grability which governs these non-perturbative ambiguities of perturbative string theory. 
This fact also makes it evidential that there is a universal strong-string-coupling dual 
theory which tightly connects various non-critical perturbative string theories in non- 
perturbative regime. 

At the first sight, the multi-cut boundary condition still seems like artificial constraints 
which specialize the solutions (in string equation) originated from the matrix models. In 
particular, the physical origin of the complementary boundary conditions is still unclear. 
In the future study, we should clarify this point and uncover the physical origin of these 
constraints. There is, however, an interesting possibility to solve this problem. It is about 
a continuum description of (/, fc)-FZZT branes. From the CFT approach, the boundary 
states are described by a superposition of (1, 1)-FZZT branes with a shift of boundary 
cosmological constants [5HI55]: 

T;{l,k)) = V|r+(- + -)27r;(l,l)) , C = v^cosh(pr). (4.1) 

/ FZZT I ^ p I FZZT 

The description in the matrix models was also found and studied in [97] . From this point 
of views, the general interpretation is that the (1, 1)-FZZT branes are fundamental and 
the others are expressed by them and that this is a reason why we have a single resolvent 
operator which describe the fundamental one. However, the above formula also suggests 
the following thing: The shifts of the momentum in Eq. (14. ip are essentially related to 
choosing different (independent) solutions of the Baker-Akhiezer system fl2.14p on which 
we impose the complementary boundary conditions. Therefore, this implies the following 
possibility: it is true that all-order perturbative amplitudes of {I, k)-FZZT branes are 
equivalent to those of the (1, 1)-FZZT branes with the relation however they are 

described as non-perturbatively different objects in string theorycj That is, the above 
relation is accidental in perturbation theory and these (/, A;)-FZZT branes have different 
realizations in the matrix models. If this consideration is true, then this would mean that 
the non-perturbative completion of string theory requires non-perturbative realization of 
all the D-brane spectrum in the theory. 

Another approach to the physical meaning would be theoretical consistency of non- 
critical M theory, i.e. the requirement that non-critical M theory can be consistently re- 
constructed. It is because we expect that the non-perturbative ambiguities should have 

After publication of this paper, we were informed that it was aheady observed that the Seiberg-Shih 
relation (|4.ip does not generally hold beyond the disk-amplitude level (i.e. even at the perturbative 
level) [TOOlfTOT] . We would like to thank Max R. Atkin for drawing our attention to the papers. 
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Non-critical M theory (r > 1) 


M = 2 QFTs (r = 1) 


perturbative string-theory sectors 


local vacua dxW{x) = 


ZZ-branes connecting string theories 


solitons connecting local vacua 


Stokes multiphers (D-instanton fugacity) 


Stokes multipliers (the number of solitons) 



Table 2: The dictionary of isomonodromy terminology between non-critical M theory and Af = 2 QFTs. 
Local vacua in J\f ~ 2 QFTs are interpreted as perturbative string-theory sectors in non-critical M theory 
(of weak coupling limit). Solitons among the vacua are understood as ZZ branes among the perturbative 
string-theory sectors. Interestingly, the Stokes multipliers in A/" = 2 QFTs are interpreted as physical 
observables. Therefore, this suggests that the D-instanton chemical potentials in string theory should 
be naturally calculable as physical observables in non-critical M theory. 

definite physical meaning in the strong-coupling dual theory, although non-perturbative 
ambiguities cannot be fixed in the perturbative frameworko An interesting clue to this 
approach is A/" = 2 QFT of Cecotti-Vafa [98j. Interestingly, there appear the k x k 
isomonodromy systems of Poincare index r = 1 and therefore one can find a translation 
into non-critical string/M theory (shown in Table [2]). As one can see, non-critical M 
theory is a simple generalization of A/" = 2 systems to the general Poincare index r. The 
critical points of r = 1 (i.e. g = 0) in the multi-cut two-matrix models are lower critical 
points than the first critical points (r = 2) which is the "pure-gravity" counterpart in 
the multi-cut two-matrix models. This situation is an analogy of the relation between 
topological series {p,q) = {^,q) in minimal string theories and the other critical points 
which start with the pure-gravity {p,q) = (2,3) points (i,e. p>2). In this sense, the re- 
construction of the multi-critical (r > 2) non-critical M theory is to find the multi-critical 
generalization of A/" = 2 QFTs. Also in this reconstruction, it is interesting to see how the 
quantum integrability of the Stokes multipliers play its rule. It is also interesting if there 
is a connection between quantum integrability of the multi-cut boundary conditions and 
the one of Nekrasov-Shatashvili |99j . 

The appearance of the quantum integrability in the Stokes multipliers of the general 
isomonodromy systems also open interesting directions. In particular, since the spectral 
parameter in non-critical string theory is understood as the spacetime coordinate [25|l36] , 
it is interesting if the new spectral parameters from the T-systems are understood as the 
third dimension of non-critical M theory. 
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A Sketch of a proof for Theorem [1] 

Here we show a brief sketch of a proof for Theorem [H The procedure is the following: 

1. Introduce a division the vectors (Eq. (13. 9 p for r = 3) and the equation Eq. (13. 5p of 
n = according to modulo r. 

2. Write down all the equation Eq. (13. 5 p of n = according to the modulo-r division. 

3. Obtain a closed set of equations, by using proper shift rules. 

4. Solve them and obtain equations with non-zero Xn'^\ With identifying the sequence 
{Un^^} with non-zero Xn'^\ these equations result in Eq. (I3.18p . The shift rules of 
the equation then generate all the equations of Eq. (13. 5p . 

Here we demonstrate this procedure in the case of {k, r) = (6m -|- 1, 3). 



Step 1 Here we introduce a modulo-r division of the vectors Eq. (13. 9 p and Eq. (13. 5p 
for n = 0. The division is given in the following way: 



4-, X, 



i 

4^, X 



X, (3m + 1)' 



X 



, (3m + 1)' 



There are only non- 
' zero components, i^. 



X, X, (3m + 4)^ 


; 


J 


k 2 ^ k X ^ k 


1, 


2, 


^ ) 




i ■ 


■ i 




i 


i 


i 


5 5 


5 


5 


k 2^ k X^ k 


1^, 


X, 


^ ) 



The number of the components is a multiple of r. 



(A.i: 



Here the components (1, 2, 3) are moved to the next to the component "fc" . The meaning 
of the arrow is 



(sym) ^(6) 



X . A . 



(A.2) 



Here "A" of Xjl means that the x^l 7^ 0. One can always prove that the contributions 
with Stokes multipliers are always accompanied with the non-zero component x^^^^ . Then 
we first divide them into a part which consists only of non-zero components and the 
other part, in which the number of components is a multiple of r. We introduce modulo- 
r division in the latter part so that the right-hand side of k should be boundary of the 
division, "■ ■ ■ , k\l, 2, 3)". In this way, the equations below become simpler. 
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Step 2 We write all the equation Eq. (13. 5p for n = 0: 

4^, X 



I. 



4^, X, X, {3m+iy 

••, X, (3m +1) 



A 



SO) 



„{6) 



(3n+l)^ ^{3n+l)'^' ^ 



II. 



III. 



1, 2, 3 
1^, X, X 



= + e^x^'l 



„{6) , /j(O) Jo) 
^2 ' 2 4 
_ ^(6) I o(0) { 



A; — 3n + 1, /c — 3r;, + 2, A; — 3n + 3 
k — 3n + 1, k — 3n + 2, k — 3n + 3 



= 1, 2, ■ ■ ■ , m — l) 



IV. 



X , X , (3m + 4) 

3m + 2, 3m + 3, 3m + 4 



A -\ 



{k-3m + 3) 



A -\ 



^ < 



(0) _ (6) .ni^^'^-^)) 
^3m+4 ^{3m+4)^ '^6m-2 "^(Sm+l) 

(0) _nb(^-^))m 

•^3m+3 ~ f^em-l ■^{3m+l)A 

(0) _ n{<^(^-^)) (0) 
■^3m+2 '^Gm (Sfh+l)^ 



k — 3m + 1 k — 3m + 2, A; — 3m + 3 

(0) 

(A.3) 



Step 3 We make a shift rule to obtain equations among 



{6n) (6n) (6n) 



III. <! J6n) _ J6{n+1)) .(6(2n-l)) (6„) ^(6(2n)) (6n) ( ^ _ . n 3) 



{6m) _ ^(6(m+l)) /i(6(2ni-l)) (gm) 



IV. <^ 



^3 ^3A I >^6m-2 •^(6m+l)A 

{6m) _ (6(2m-l)) (6m) 

•^2 — '^6m-l -^(Gm+l)^ 

(6m) _ ^(6(2m-l)) (6m) 

•^1 "^em •^(6m+l)'^ 



(A.4) 



Step 4 By solving these equations, we obtain equations for non-zero components of 

in— I 



U - "t.tA + ^(in~2 -^(l+en)^ 



^ (6(2n)) (6„) 

C'6„+i -^(4+6,1)^' 



71=1 



71=0 



m / \ 
n _ \^fll^(2n-l)j (6„) 
^ ~ / .^&n-l ■^(6n+l)A 



m—1 



X fr- , T \ A ~t~ 



(6(27i)) (6n) 
6n+2 ^(6n+4)A' 



n=l 



n=0 
m—l 



m / \ m—i / \ 

(0) _ l6{2n-l) j (6„) a I6(2n) j (671) 

1^ ~ Z^^67i ■^{l+6n)A / . ^^671+3 •^(4+6n)A- 



n=l 



n=0 



By introducing the following sequence {vn^^ 



,(0) ^ ^{6n) 



— ^(371+1)^ I ^(3n+l)^ 



,(6(n-l)) 



X 



(6(71— m)) 
(3n+l)^ 



(A.5) 



(A.6) 
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we obtain the equations of the form: 



2m 



^1 



^2 



(0) 
y2fh 



3) (0) _ p 

m+1 ^ '^3(n-l)+l ~ ^' 

n=l 

2m 

(0) _ g 



n=l 



2m 



-1/6 



(0) 



V- ^(6(n-l)) (0) _ Q 



(A.7) 



n=l 



By taking into account Eq. (13.171) . we obtain Eq. fl3.18p with s = 0. Therefore, all the 
other equations in Eq. fl3.18p are obtained by the 2s shift rule. Note that yn^^ is defined 
as 



yn 



(2s) 



X 



{6n+2s) 



(A.8) 



In this sense, we use the notation of k = rm + / (1 < / < r) in the following discussions. 
Finally we mention the extension of this procedure to all the {k, r) cases. We note is that 
we have used a specialty of r = 3, i.e. mi = 1 (See Eq. (13.11) ). Because of this specialty, 
the division procedure Eq. ( lA.ip is simpler and the proof is not so complicatedo In the 
general cases, however, we can re-organize the components of the vector X^"^ as 



1 ' "^l+rmi ' '''l+2rmi ' 



(0) 



(A.9) 



In this way, we can perform the same procedure in general r and can show Eq. (13.181) for 
general {k,r). 
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